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Lieber Derive Anwender 


Dear Derive User, 


Die erfreuliche Nachricht zuerst: die DUG wachst 
weiter. Die Mitgliederzahl ist inzwischen auf iiber 
180 angewachsen. Voraussichtlich wird sich die 
DUG noch heuer weltweit etablieren. Die Produk- 
tentwikkler von DERIVE sind mit einem dies- 
beziiglichen Vorschlag an uns herangetreten. Wir 
sehen mit Interesse der Entwicklung entgegen. 

AuBerdem mochte ich Ihr Augenmerk auf die 
Moglichkeit lenken, DERIVE Consultant von SWHE 
zu werden. Ein Anmeldeformular liegt dieser Auflage 
bei. 

Viele interessante Zuschriften haben den Her- 
ausgeber wahrend des Sommers erreicht. Drei davon 
mochte ich besonders herausgreifen: Mr. Setif aus 
Annemasse hat mir insgesamt iiber 100 Manuskript- 
seiten iiberlassen. Sie sind voll mit Ideen, was man 
mit DERIVE alles behandeln konnte: Ketten- 

bruchentwicklungen, Primzahlprobleme, Folgen und 
Reihen, Darstellung von Polygonen, Fraktale und 
noch viel mehr. Der D-N-L hat nun deshalb fur ei- 
nige Zeit eine eigene Kolumne "Mr. Setif s Treasure 
Box" eingerichtet, um seine Ergebnisse und Anre- 
gungen zu publizieren. Zwei Beispiele konnen Sie 
mit 

SYRACUS.MTH und mit einer Auswahl von Frakta- 
len auf Seite 7, bzw. auf der inneren Umschlagseite 
finden. Der zweite umfangreiche Beitrag, mit dessen 
Ausgabe wir im nachsten D-N-L beginnen werden, 
stammt von Prof. Hodgkinson, Liverpool, und be- 
schaftigt sich mit dem Einsatz von DERIVE beim Er- 
arbeiten des Folgen- und Reihenbegriffs. Dieser Arti- 
kel scheint mir deshalb besonders bemerkenswert, 
weil er viele methodische und didaktische Elemente 
beeinhaltet. Eine Zuschrift der Herren Appel und 
Brand zeigt neue Moglichkeiten, allerdings in einem 
geometrischen Zusammenhang auf. 

Die nachste Ausgabe wird wesentlich vielfalti- 
ger sein als diese, deren Einheitlichkeit dadurch be- 
dingt ist, daB ich den Beitrag iiber die Differential- 
gleichungen abschlieBen wollte. Ich wiinsche Ihnen 
trotzdem viel Vergnugen mit D-N-L#3 und freue 
mich iiber jede Zuschrift, Anregung und/oder Kritik. 

AuBerdem werden Sie vielleicht im D-N-L#4 
eine Uberraschung vorfmden! 


Let me first of all bring you good news: DUG keeps 
on increasing and has meanwhile been joined by over 
180 members. Most probably DUG will be making 
its way world-wide before the end of the year. As a 
matter of fact, the developers of DERIVE have ap- 
proached us with a proposal relating to such an ex- 
pansion. Thus we are eagerly awaiting further deve- 
lopments in this matter. 

Furthermore I would like to draw your attention 
to the possibility of becoming a DERIVE Consultant 
of SWHE. An application form has been enclosed in 
this issue. 

Many interesting letters have reached the Editor 
during the summer. I take the liberty of pointing out 
three of them which, in my opinion, deserve particu- 
lar attention: Mr. Setif from Annemasse has let me 
have a script totalling over 100 pages which are full 
of ideas of what fields DERIVE may be applied to: 
continued fractions, prime number problems, sequen- 
ces and series, polygons, fractals and many more. For 
this reason D-N-L has installed an individual column, 
"Mr. Setif s Treasure Box", in which his results and 
suggestions will be published for some time. Two 
examples concerning SYRACUS.MTH and a choice 
of fractals con be found on page 7 and the inside 
cover, respectively. The second, fairly detailed con- 
tribution, which is to be published in the next issue, 
has been handed in by Prof. Hodgkinson, Liverpool, 
and deals with the application of DERIVE as being 
instrumental in establishing the concept of sequences 
and series. This article seems particularly remarkable 
to me as it comprises a good many methodical an 
didactic elements. A letter by Mr. Appel and Mr. 
Brand points out new possibilities which, however, 
have to be considered in a geometric context. 

Next issue will be considerably more varied 
than this one, its basic uniformity being mainly due to 
the fact that I wanted to complete the contribution on 
differential equations. Despite this, I hope that you 
will enjoy D-N-L#3 and would like to invite you to 
send in contributions - ideas or critical comments are 
equally welcome. 

By the way, in D-N-L#4 a surprise may be wai- 
ting for you! 
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The Derive-News-Letter is the 
Bulletin of the Derive-User- 
Group. It is published at least 
three times a year with a con- 
tent of 30 pages minimum. The 
goals of the D-N-L is to enable 
the exchange of experiencies 
made with Derive as well as to 
create a group to discuss the 
possibilities of new methodic 
an didactic manners in teaching 
Mathematics . 

Subscription of the D-N-L is 
restricted to Members of the 
Derive-User-Group . Membership- 
form is enclosed with the first 
issue you get. 

Editor : 

Mag. Josef Bohm 
A-3042 Wurmla 
D ’ Lust 1 
Austria 


Contributions : 

Please send all contributions 
to the above adress. Non-Eng- 
lish speakers are encouraged to 
write their contributions in 
English to confirm the interna- 
tional touch of the D-N-L . At 
the other hand non-English ar- 
ticles will be warmly welcomed, 
too. Your contributions will be 
edited but not refereed. By 
submitting articles the author 
gives his consent for reprin- 
ting in the D-N-L. The more 
contributions you will send to 
the Editor the more livelier 
and richer in content the De- 
rive-News-Letter will be. 


Preview (Contributions for the next issues): 

With ITERATES to Chaos (Part 2) 

How to write one's own DEMO-file 
Cyclomania (Hypotrochoids and Epitrochoids) 
Joined Venture Derive and Cabri 
Laplace T ransforms 

A Module on Sequences and Series with Derive 


(will be published December 1991) 


Impressum : 

Medieninhaber : DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA; 
Richtung: Fachzeitschrif t ; 

Herausgeber: Mag. Josef Bohm, Adresse wie oben; 

Herstellung: eigene Vervielf altigung 
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Derive — User - Forum 


Dr. Felix Schumm, Stuttgart 

Zum Zeichnen von Strecken zwischen zwei Punkten xl und x2 habe ich mir in 
mein USER-File die Funktion aufgenommen, bei der nach SIMPLIFY-PLOT die 
vorgeschlagenen Grenzen des Parameters -n bis +tu nur bestatigt werden mus- 
sen . 

LINE(xl,x2) := xl + (t/n + 1) (x2 - xl) 

Aufruf z.B. mit LINE([1,2], [ — 4^3] ) 

(Felix Schumm presents a function to plot a segment between two points, which needs only pressing 
the ENTER-key in the 2D-Plot Window to accept the default values for parameter t.) 


H. Wunderling, Berlin 

. . . Mit Bewuhtsein haben Sie wohl gleich ins erste Heft den Beitrag von 
Herrn Schmidt (Kdln) getan, denn Sie 

"... erwarten auherdem den Widerspruch, der mir genehm." 

Deshalb widerspreche ich als Lehrer, der seit 2 Jahren DERIVE im Unterricht 
einsetzt. Ich bin uberzeugt, daft Werkzeuge wie DERIVE den Mathematikunter- 
richt erheblich beeinflussen werden und nicht ausgesperrt werden durfen. 

Hier aber erst einmal 3 Dinge, die ich Sie bitte, an die Entwickler von 
DERIVE wei terzugeben: 

- Der Zugang zu den einzelnen Menupunkten sollte vom Lehrer gesperrt werden 
konnen . (Ein Konf igurationsprogramm, das dem Schuler unzuganglich ist, muh- 
te dem Lehrer erlauben, genau die Teile aktiv zu halten, die er fur seine 
Schuler auswahlt. Z.B. ist der soLve-Befehl didaktisches Gift, wenn Umfor- 
mungen mit DERIVE geubt werden sollen; Oder EIGENVALUES, wenn ....) 

- Eine Schullizenz durfte preislich nicht wesentlich iiber einer Einzelli- 
zenz liegen. Es ist abzusehen, daft Computer-Arbeitsplatze fur Schuler vor- 
handen sein werden und etwas spater auch flexiblere Notebook-PCs . 

- Bei der Behandlung des zentralen Grenzwertsatzes der Stochastik habe ich 
bemerkt, daft die neue Funktion ITERATES auch syntaktisch falsch gesetzte 
Klammern zu verarbeiten beginnt; der Editor fangt die Fehler nicht ab ! 

z.B. statt: 

ITERATES (VECTOR (IF (j=l OR i=23 , 0 , 1/3 (ELEMENT (v, i-1 ) +ELEMENT (v, i ) + 
ELEMENT (v, i+1 ) ) ) , i, 1 , 23) , v, vl , 20 ) 

wird auch akzeptiert: 

ITERATES (VECTOR (IF (j=l OR i=23 , 0 , 1/3 (ELEMENT (v, i-1 ) +ELEMENT (v, i ) + 
ELEMENT (v, i+1) ) , i , 1 , 2 3 ) ) , v, vl , 2 0 ) 

Oder auch: 

ITERATES (VECTOR (IF (j=l OR i=23, 0, 1/3 (ELEMENT (v, i-1) +ELEMENT (v, i) + 
ELEMENT (v, i+1) ) , i , 1 , 2 3 ) , v, vl , 2 0 ) ) 

(Mr Wunderling demonstrates that Derive accepts parentheses which are set syntactically incorrect. In 
my opinion Derive internally counts the number of opening and closing parentheses and does not 
check, if their position make sense or not.) 
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Dipl.-Ing. K. Schmidt, Koln 

Die Reaktionen auf meine Zitate sind genau so erfolgt, wie ich es erwartet 
habe . Nach dem Auto ist der Computer die zweite "Heilige Kuh" der Menschen. 
Es ist schade, daft Sie nicht gebracht haben, was auf der Tagung uber "Op- 
portunities and Risks of Artificial Intelligence" in Hamburg 1989 festge- 
stellt wurde . Durch die Kunstliche Intelligenz werden vor allem zwei Ent- 
wicklungen verscharft: das Verschwinden der Verantwortung und die schlei- 
chende Erosion von Qualif ikation und Wissen. Einmal soil also der Computer 
alle Entscheidungen treffen und damit die Verantwortung iibernehmen, weil 
die meisten Menschen nicht denken wollen, zum Teil auch gar nicht konnen . 
Zum anderen sagen sie sich: warum soil ich mir noch Miihe geben, der Compu- 
ter macht doch alles. Sie spielen mit dem Rechner, lernen aber nicht, wie 
man es macht. Es gibt sogar Inf ormatiker , die diese Richtung propagieren. 
Schon wird in Fachzeitschrif ten vor bestimmten Expertensystemen gewarnt, 
weil sie gedankenlos genutzt werden. Ich war 19 Jahre "Knowledge Engineer". 
Man glaubt gar nicht, was fur verruckte Ideen ich manchmal programmieren 
Oder verbessern muftte, wobei es Arger mit Akademikern aller Grade gab. 


Andererseits darf man nicht vergessen, daft die Derive-User-Group ein ex- 
klusiver Klub von hochgradigen Experten ist, von denen vielleicht mancher 
nicht weift, wie es in der Praxis zugeht. Gute mathematische Kenntnisse ha- 
ben leider nur sehr, sehr wenige Menschen, die dazu die entsprechende Ver- 
anlagung haben mussen. Fur den hohen Stand von Technik und Naturwissen- 
schaften sind es viel zu wenige. Z.B. brechen bei den Studenten der Infor- 
matik (gehort ja zur Praktischen Mathematik) ca 70% (1989 waren es in 

Frnnkfurt sogar 88%) das Studium ab . 


Robert Setif, Annemasse, Frankreich 

1. Voici quelques fichiers (files) *.MTH. Beaucoup de mes fonctions sont 
lentes et lourdes inutilement, mais elles sont plutot destinees a des 
apprentis en programmation . En effet les exemples des fichiers *.MTH de 
DERIVE ne sont pas tou jours faciles a comprendre avec un niveau mathema- 
tique insuffisant. 

2. Je commence un gros fichier sur les fractions continues pour des appro- 
ximation de V 2, V 3, .... Je sais bien que DERIVE peut les calculer di- 

rectement, mais c'est pour le plaisir de la programmation. Ces fractions 
continues vous interessent, elles? 

3. DERIVE 2 . xx est vraiment WONDERFUL avec sa nouvelle capacite de program- 
mer. Cependant il me semble qu'il manque un equivalent de WHILE conditi- 
on ou UNTIL condition, ce qui permettrait de stopper une suite d’ itera- 
tions dont on ne connait pas le nombre a 1 ' avance (voir par exemple dans 
SYRACUS.MTH) . II manque aussi une fonction qui calculerait des nombres 
aleatoires (RND) e N ou e Q ou e R. II faudrait de meme une fonction 
MOD dans N ecrite en langage machine plus rapidement que celle de 
MISC.MTH. Ce dont bien etre possible puisque NEXT_PRIME existe. 

Pourriez vous transmettre ces 3 revendications "syndicales" a Hawaii? 

J' imagine que vous en avez deja d'autres a envoyer? 

4 . Peut-etre pourriez-vous aussi publier dans votre bulletin europeen de 
DERIVE des articles concernant MuMATH? 

A propos de MuMATH, je n' ai pas re^u de Newsletter MuMath depuis novem- 
bre 88. 






D-N-L #3 


Derive 


User 


Forum 


P 5 


D-N-L : To the first part of this letter see the Letter of the Editor! 
I've sent Mr Setif's wishes to Hawaii and here is Mr. A. Rich's answer: 
good suggestions. They will be added to the Version 3 wish list." 
Maybe there are some DERIVE-Users who would like to correspond with 
Mr Setif about MuMath. I would be glad to knot the wires. 


Std. W. Kramp, Berlin 


Im VGA-Graf ikmodus wird bei mir 
die aktuelle Funktion falsch ab- 
gebildet (siehe Bild) . Wie lafit 
sich der inverse Ausdruck der ak- 
tuellen Funktion abstellen Oder 
der Fehler beheben? 



COffilAND: Build Calculus Declare Expand Facto 

Options Plot Quit Benov e Sinplifg Transfer : 


D-N-L : Another answer from A. Rich, concerning one question of Mr Lauer in 
D-N-L#2, page 6: "memory full" in GCD (123456, 987 6) : 

"due to stack overflow caused by an excessive number of recursive function 
calls" . 



Here is an example of a "needle 
graphic", produced by connecting 
the points of three curves. 

(Plot Options for the points: 
Discrete Large) 




UECTOR 


2 COS(t) - 0.5 
0.5 SIN(t) 
COS(t + 1) 


- SIN(t) 
COS(t) + 1.5 
0.5 SIN(2 t) 


t , 0, 2 n 


20 
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/oft Warehouse: 

SCHLOSS HAGENBERG CUfOOG ' 
A-4232 HAGENBERG, AUSTRIA ■ 

phone (0)7236/32 97-81, fax (0)7236/32 97-30 


WOULD YOU BE INTERESTED IN BECOMING A DERIVE CONSULTANT? 


Soft Warehouse Europe (SWHE) , the European representative of the de- 
veloper of DERIVE, is looking for knowledgeable DERIVE users who ha- 
ve experience in giving lectures, seminars, or workshops on DERIVE 
and/or in writing articles about the system. If you meet these re- 
quirements and if you would like to offer your expertise to SWHE, 
authorized DERIVE dealers, conference organizers, publishers, and 
others, we would welcome your application to become a DERIVE consul- 
tant . 


We offer to a DERIVE Consultant: 

- listing in the "European DERIVE Consultant Directory" - this 

document will be regularly distributed to authorized DERIVE 
dealers and will be sent to people or institutions requesting 
it (conference organizers, publishers, etc.), 

- a free upgrade service for the DERIVE software and manual - 

including all non-announced, minor improvements on the system, 

- one year free membership to the DERIVE User Group - including the 

receipt of the DERIVE News Letter, 

- regular information about DERIVE-related publications and presen- 

tations . 


We accept from a DERIVE Consultant: 

- active promotion of DERIVE through lectures, seminars, workshops, 

articles or books, 

- the readiness to cooperate with authorized DERIVE dealers, confe- 

rence organizers, publishers, etc. (honorarium to be agreed 
individually between the partners) , 

- the readiness to speak to potential new users who are seeking an 

independent opinion about DERIVE (authorized DERIVE dealers 
may give prospective customers/purchasers names and addresses 
of their nearest DERIVE consultants) . 

We further expect the applicant to be a registered DERIVE user. 


How to become a DERIVE Consultant: 

Return the enclosed DERIVE Consultant Application Form together with 
all requested documents. 
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A Famous Sequence 
(Working with ITERATES) 


In original DNL#3 the MOD-function had to be preloaded from MISC.MTH (together with FLOOR and 
many others which became later part of the Derive core. This is the MOD-function from those times: 


MOD (a, b) 


: = lim 
x^>a+ 


b • AT AN 


COT 


ttl: 

#2 : 
#3: 


sy(n) := 

If MOD(n, 2) 
n/2 

3 n + 1 


sy(B) 

sy(8) 


16 

4 


The element following number n: 


#4: 

#5: 
#6 : 


nsy(n) := 

If n = 1 
0 

If sy(n) = 1 
1 

1 + nsy(sy(n)) 
nsy(15) = 17 

syr(n) := ITERATES([1 + u 


The number of iterations which are necessa- 
ry to come to an end of the sequence initiali- 
zed by n. 


sy(u >1,. v, 

2 I 


[0, n], nsy(n)) 


#7: syr(3) 


0 

1 

2 

3 

4 

5 

6 
7 


3 

10 

5 

16 

8 

4 
2 
1 


The 3n+1 Sequence starts with a number n e N. 
If n is even you may halve it, otherwise multiply by 
three and add 1 . Perform the same with this newly 
created element and go on. The sequence leads to 
a power of 2. 


The complete sequence: 

#9: syra(n) := ITERATES(sy(a), a, n, nsy(n)) 


ttl0: 

syra{5) = | 

[5, 

16, 

8, 

mi: 

syra{15) = 

[15 

, 46 

, 2 

#12: 

DIM(syra(27)> 

= 112 


maxsi|{n) := ITERATES< 

f 1 

#13: 




l 



0 

34 

1 



1 

17 

34 



2 

52 

34 



3 

26 

52 



4 

13 

52 



5 

40 

52 



G 

20 

52 

#14: 

maxs</{34) = 






7 

10 

52 


sy(u }, MAX(u , u >1, u, [0, n, 1], nsy(n>) 
2 2 3 J 


The matrix shows the maximum element after k 
steps (1st col) in the last column. 

In the next issue I will show some improvements 
Mr. Setif has asked for! 
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With ITERATES to CHAOS 

Dr. Felix Schumm, Stuttgart 

Inhalt: 

1 Einleitung 

2 Iteration mit einer Variablen 

2. 1 Allgemeine Betrachtungen 

2.2 Umsetzung in DERIVE 

2.3 Beispiel (logistische Gleichung) 

2.4 Graphische Darstellung mit Hilfe der Befehle VECTOR und ELEMENT 

3 Iteration mit zwei miteinander verkoppelten Variablen 

3 . 1 Allgemeine Betrachtungen 

3.2 Umsetzung in DERIVE 

3.3 Beispiel (Volterrazyklus) 

3.4 Graphische Darstellung 


1. Einleitung 

Im folgenden Beitrag soil gezeigt werden, wie man die Befehle ITERATES, VECTOR und ELEMENT 
giinstig zur Losung von Iterationsproblemen einsetzen kann, und wie die erhaltenen Ergebnisse mit 
DERIVE graphisch dargestellt werden konnen. Es wird versucht, von einfachen Beispielen ausgehend 
zu komplizierteren Befehlskombinationen langsam fortzuschreiten. 


The following contribution shall demonstrate how to apply ITERATES, VECTOR and ELEMENT for 
solving iteration problems. We will try to proceed from easy examples to more complex combinations 
of commands. ELEMENT is still available in recent Derive versions but can be replaced by the sub- 
script command sub (i). 

2. Iteration mit einer Variablen 


2.1 Allgemeine Betrachtungen 


Einer Zahl x a i t werde mittels eines Funktionsterms/jxait) ein neuer Wert x neu zugeordnet. 

/■ 

X neu=f( X alt\ bzW ' X neu X alt 

AnschlieBend wird die Zahl x neu mit dem gleichen Funktionsterm erneut abgebildet, bzw. als ein neues 
Xait interpretiert. Fiihrt man diesen Prozess zz-mal durch, so spricht man von zz-facher Iteration. 


Fur n = 3 gilt: 


/: /: /: 

X neu = f(f(f( X alM b ZW. X mu <- X <- X <- X aU 


Programmiert wird eine Iteration durch eine zz-fache Schleife an deren Ende der Inhalt der Variablen 
-^neu jeweils der Variablen x a i t zugewiesen wird. 


LET X a it • -^neu 


A X alt) 

A 


Schleife n - mal durchlaufen. 

n-fold Loop 
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In BASIC lautet ein entsprechendes Programm etwa so: 

xalt = startwert 
FOR i = 0 TO n 

PRINT " Iterationsstuf e : "i"; x-Wert: n xalt 
xneu = f (xalt) 
xalt = xneu 
NEXT i 

Man benotigt also einen Funktionsterm/(x), einen Startwert xo fur die Variable x und die 
Iterationstiefe n. 


2.2 Umsetzung in DERIVE 

Die Iterationsschleife ist in DERIVE kompakt programmierbar mit dem Befehl ITERATES. 


Es sei: 


Nach Eingabe von 


X 

die Variable, die iteriert wird, 

f 

irgend ein Term in x, 

xO 

der Startwert von x und 

n 

die Iterationstiefe. 

von 

F(X) := 


ITERATES (F (X) ,x,xO,n) 


liefert SIMPLIFY (oder APPROXIMATE) die Ergebnisse der einzelnen Iterationen in Form eines 
Vektors mit folgender Bedeutung der Komponenten: 


xO 

1 

1 

Startwert ( Initial value ) - keine Iteration (no iteration ) 

f (xO) 

1 

Xneu (fur 1-mal iteriert) 

f (f (xO) 

1 

1 

Xneu (fur 2-mal iteriert) 

f (. . . (f (xO) . 

1 

• .) 1 

Xneu (fur n-mal iteriert) 


J 



2.3 Beispiel 

Fiir die Wachstumsrate einer Population gilt unter bestimmten Voraussetzungen folgende 

Differentialgleichung : ( Differential equation for modeling the growth rate of a population ) 

dx 7 2 

— = a x-b x 
dt 

t : Zeit - time 

x: momentane Anzahl von Individuen (bzw. momentane Dichte von Individuen je 

Flachen- oder Volumseinheit) 

instantaneous number of idividuals (density of individuals per unit of area or vo- 
lume) 

a x Vermehrungsrate proportional zur Anzahl der schon vorhandenen Individuen 
growth rate is proportional to number of existing individuals 
b x 2 : Abnahme durch gegenseitige Behinderung bei groBer Individuenanzahl. 

decrease because of a negative mutual influence fora large number of indivi- 
duals 
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Die Differentialgleichung laBt sich leicht integrieren und ergibt eine Losung, die fur t — > +co immer 
die Asymptote x = alb hat. (DNL #2, S. 25). LaBt man die Zeit (und damit auchx!) nicht kontinuierlich 
anwachsen (biologisch bedeutet dies, daB nur ganzzahlige Nachkommen moglich sind), so andem sich 
die Losungsmoglichkeiten drastisch und man erhalt ein Iterationsproblem. Aus der Differentialglei- 
chung wird zunachst die zugehorige Differenzengleichung: 


The DE can easily be integrated and gives a solution, which for t — > +oo has always y = a/b as asym- 
ptote. (DNL#2, p 25). Taking discrete steps for t and x - instead of contineous ones - which has the 
biologic interpretation that only an integer number of succeeders is possible the variety of solutions 
changes dramatically and we are facing an iteration problem. First of all we make a "difference equa- 
tion" from the given differential equation: 


— = a x-b x 2 wobei Dt ein endliches Zeitintervall bedeutet. 

Dt 

Dx = [a x-b x 2 ) • Dt 

X neu= X a„ +Dx 

X neu= X al, + ( a - X alt +b - X au)- Dt 
X neu=Q- + a - Dt )- X al,- b - Dt - X aU 

Der Funktionsterm, den wir untersuchen wollen, lautet also: 

f (x) := (l+a*Dt) *x-b*Dt*x A 2 

Zur Vereinfachung wahlen wir das Zeitintervall Dt = 1 / (b-a ) und erhalten: 


2 

f(x) := (1 + aDt)x - b Dt x 

f 1 1 12 

f(x) := 1 + a x - b x 

[ b - a J b-a 

2 

b -x b -x 


a - b b-a 


— >f (x) : =r ( x-x A 2 ) mit 


b 

r = 

b-a 


Zeichnet man sich die Zuordnung x neu =f{x a i t ) auf, so ergibt dies eine abwarts gerichtete Parabel mit 
den Nullstellen (0|0) und (1|0) und dem Scheitel (0.5|r/4). Man erkennt, da!3 nur Werte 0 < x a i t < 1 fur 
den Start sinnvoll sind, da sonst schon im ersten Iterationsschritt die x neu - Werte negativ oder 0 wiirden, 
was fur unser Problem keinen Sinn macht. 

Man kann die Frage stellen, ob es Startwerte gibt, bei denen die Individuenanzahl vollig unverandert 
bleibt. Fur diese nicht notwendig stabilen Gleichgewichtswerte oder Fixpunkte der Abbildung/{x) gilt: 

The sketch ofx neu = f(x a / t ) shows a downward directed parabola with zeros (0|0) and (1 10) and vertex 
(0.5|r/4). It is obvious that only values 0 < x ait < 1 make sense, because otherwise the first iteration 
would result in 0 or in negative values. One can ask now it there are initial values which leave the 
number of individuals unchanged. So we are looking for the fixponts: 


Xfix =A X fix) und damit x flx = r x fix ( 1 - x flx ) 


(und den Trivialfall x fix = 0). 


und man erhalt sofort nach Auflosen dieser Gleichung: x 


r 
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Naheliegenderweise interessieren uns nur positive Xfi x -Werte, d.h. wir werden nur r-Werte mit \r\ > 1 
untersuchen. 


Nun zur Programmierung in DERIVE. Als Voreinstellung wahlen wir: 

Option Notation (decimal, digits 3) 

Option Precision (approximate) und geben ein: 


#1: 

Precision := Approximate 

#2: 

Notat ionDigits := 3 

#3: 

2 

f := r (x - x ) 

#4: 

r - 1 

x_fix := 

r 

#5: 

ITERATES(f, x, x0, n> 

#6: 

[n := 13, r := 2] 

#7: 

x_fix = 0.5 


Wir berechnen den Gleichgewichtswert fur r = 2 und erhalten 0.5. 

Als Startwert definieren wir willkurlich 0.1, markieren die Zeile #5 und fiihren Simplify durch. Dann 
erhalt man die einzelnen Stufen der Iteration. 


#8: x0 := 0.1 

#9: [0.1, 0.18, 0.295, 0.416, 0.485, 0.499, 0.499, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5] 

#10: x0 := 0.9 

#11: [0.9, 0.18, 0.295, 0.416, 0.485, 0.499, 0.499, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5] 

Wir konnen aber auch einen anderen Startwert wahlen (z.B. xO = 0.9) und fiihren anschlieBend wieder 
Zeile #5 aus: Gleichgiiltig, welchen Wert man fur xO (0 < xO < 1) ansetzt, immer lauft die Iteration auf 
den Gleichgewichtswert Xf lx = 0.5 hin. 

We can choose several different initial values xO, approximating or simplifying expression #5 always 
converges to x fix = 0.5. 


Dies ist jedoch nicht fur jedes \r\ > 1 der Fall! Wahlen wir z.B. r = 3.2, so ergibt sich x flx = 0.687. Die- 
ser Fixwert wird jedoch nicht angenommen, vielmehr schwanken die Ergebnisse, wenn man lange 
genug iteriert, unabhangig vom Anfangswert xO zwischen den beiden Zahlen (= Attraktoren) 0.799 
und 0.513. Der Zuwachs, bzw. die Abnahme von x schieBen jeweils gewissermaBen iiber den Gleich- 
gewichtszustand x flx hinaus. Die entsprechenden Eingaben in DERIVE sind (fur die erste Iteration gilt 
noch immer xO = 0.9): 


#12: 

r := 3.2 















#13: 

[0.9, 0.288 

!, 0. 

656, 

0.721, 

. 0.642, 0.735, 0.623 

, 0.751, 

F 0.597, 

. 0.769, 

0 

.567, 

0. 

.785, 


0.539, 0. 

795] 














#16: 

ITERATES(f 

.. x. 

0.1 

i- n) 

= 

[0.1, 

0.288 

, 0.656, 

0.721, 

0.642, 

0.735, 

0. 

623, 

0. 

751, 


0.597, 0 

.769 

, 0. 

567, 

0. 

.785, 

0.539, 

0.795] 








#17: 

ITERATES(f 

.. x. 

0.5 

- n) 

= 

[0.5, 

0.8, 

0.512, 0. 

.799, 0 

.512, 0. 

.799, 0. 

.513, 0. 

.799, 


0.513, 0 

.799 

, 0. 

513, 

0. 

.799, 

0.513, 

0.799] 








#18: 

ITERATES(f 


0.8 

* n) 

= 

[0.8, 

0.512 

, 0.799, 

0.512, 

0.799, 

0.513, 

0. 

799, 

0. 

513, 


0.799, 0.513, 0.799, 0.513, 0.799, 0.513] 
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Was geschieht fur xO = x flx ? 

#14: 0.687 

f 3.2-1 ] 

#15: ITERATES f, x, , n = [0.687, 0.687, 0.687, 0.687, 0.687, 0.687, 0.687, 

l 3.2 J 

0.687, 0.687, 0.687, 0.687, 0.687, 0.687, 0.687] 


Die Begriindung, warum fur r > 3 der Gleichgewichtszustand x flx = (r-l)/r unstabil wird, soil hier 
nicht weiter verfolgt werden, da man diese in vielen Artikeln zum Thema CHAOS unter dem Schlag- 
wort "Feigenbaum" (oder "Bifurkation") nachlesen kann. Erwahnt sei nur, da6 fur 3 < r < 3.4495 zwei 
Attraktoren, fur 3.4495 < r < 3.75 vier und fur r > 3.75 acht und mehr Attraktoren auftreten. 

We will not explain here why for r> 3 the equilibrium x fix = (r-1)/r becomes unstable, because there is 
so much reference literature about CHAOS ("bifurcation", "Feigenbaum"). It shall be mentioned that 
for 3 < r< 3.4495 two attractors are appearing, for 3.4495 < r 3.75 there are four and for 
r> 3.75 there are eight or more attractors. 


2.4 Graphische Darstellung 

Gibt man eine Matrix, wie z.B. 


2 3 
4 5 

3 7 


in DERIVE durch: [[2,3], [4,5], [3,7]] (oder spater durch [2, 3; 4, 5; 3, 7]) ein, so liefert 
der 2D-Plot ein Schaubild, in welchem die Punkte (2|3),(4|5) und (6|7) eingetragen sind. (Im Options 
State muss man Rectangular und Discrete wahlen - Derive for DOS!). Um unsere Iterationen gra- 
phisch darzustellen, sollten wir aus dem Ergebnisvektor 

v := [xO, f (xO) , f (f (xO) ) , ] 


eine Matrix "Bild" mit folgender Form herstellen: 



0 


xO 


0 

ELEMENT 

(v, 

1) 1 


1 


f (xO) 


1 

ELEMENT 

(v, 

2) 


2 

f i 

( f (xO ) ) 


2 

ELEMENT 

(v, 

3) 

Bild : = 

n 

f i 

(f(f(. .). 

Bild : = 

n 

ELEMENT 

(v, 

n+1) . 


0 

1 


Bild : = 


2 


v 

1 

V 

2 

v 

3 


n v 

n+1 


Mit dem ELEMENT-Befehl konnen wir die Matrix Bild aufbauen und mit Hilfe von VECTOR laBt 
sich diese Matrix iiber eine Schleife erzeugen: 


BILD : = VECTOR ( [i, ELEMENT (v, i+1) ] , i, 0,n] 
oder mit Derive 5 

BILD : =VECTOR ( [i, v SUB ( i+1 ) ] , i , 0 , n] 


Markiert man diese Zeile und fuhrt Simplify aus, so erhalt man genau solche Matrizen, die unsere 
Iterationsergebnisse geeignet fiir das PLOT-Fenster aufbereiten. Mit Hilfe der Tasten F3 und F4, so- 
wie mittels Substitute lassen sich die Ergebnisse bequem zu einer Funktion "BILD” zusammenfassen. 
Der unnotige Rest wird geloscht. Wir konnen z.B. ausschlieGlich folgende Zeilen beibehalten und in 
einem File speichern. 
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ttl: [r :=, x0 : =] 

#2: F(x, r) := r x (1 - x) 

#3: ergeb_uektor(r, x0, n) := ITERATES(F(x, r), x, x0, n) 

BILD(r, x0, n) := UECTORCfi, (ergeb_uektor(r, x0, n)> 1, i, 0, n) 

#4: L i + lj 


Geeignete Werte fur (r,xO,n) sind z.B. (1.2,0.001,100), (1.2,0.4,100), (3.2,0.687,100), (3.5,0.714,100) 
oder (3.6,3.722,100). So liefert die Eingabe: 


#5: BILD(3.2, 0.687, 100) 

nach Simplify die folgende Bildmatrix: 

0 0.687 

1 0.688 

2 0.686 
3 0.688 


98 0.513 

99 0.799 

100 0.513 


#7: BILD(1 .2, 0.001, 100) 


#8: BILDC1.2, 0.4, 100) 


Mit Plot kann das entsprechende ^-x-Diagramm 
ausgegeben werden. 





#9: 


BILD(3 .2, 0.687, 100) 
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#10: BILD<3.5, 0.714, 100) 


#11: BILD(3 .6, 0.722, 100) 




wird fortgesetzt (to be continued). 



A Family of Trochoids 
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Solving Odes Using 

Derive ' s 0DE1.MTH, Part 2 

Josef Bohm, Wiirmla 


I tried to extend the file ODE1.MTH in orderto be able to write general solutions with a parameter c 
instead of (xOj^O) used in ODE1.MTH. You will find a function to express the integrating factor and a 
function to solve the LAGRANGE-DE. The expressions #1 - #4 have been used in Part 1 and will be 
used now, too. I collect all the new functions and hints in the file ODE1EXT.MTH. Having loaded 
ODE1, 1 merge (load) ODE1EXT and start working. 


(In later MATH-files ODE1.MTH general solutions using a constant c have been added.) 


ODE1 EXT.MTH 


#1 : GEN_SOL_SEP(p, q, x, y) 


T dy - J p dx = c 

J q 


#2: SPEC_SOL_SEP(p, q, x0, y0, x, y) : = 


y 

I - 

J q 


dy 


x0 


p dx 


#3: 


DIR(r, x, y, x0, y0) := lin lin IF 

y-*y0 x-*x0 


y0 

1 


= 0, [x, y + t]. 


sT(l + r ) 


t r 


2 

sT(l + r ) J 

#4: FELD(r, x, y, xl, xr, xs, yd, yu, ys) := UECTOR(UECTOR{DIR(r, x, y, x0, y0), x0, 

xl, xr, xs), y0, yd, yu, ys) 


#5: 

LINEARl(p, q, x, y, x0, y0) 


y 


x 

P ^INT{p, x, x0, x) 
y0 + J q ■ e dx 

x0 


^INT(p, x, x0, x) 
e 


ttG : 

U?z 

#8: 

#9: 
# 10 : 
# 11 : 

#12: 


GEN_LINl(p, q, x, y) 


y = 



p dx 


dx + c 


J 


P 


dx 


GEN_BERN(p, q, k, x, y) := GEN_LIN1{{1 - k) p, 

GEN_S OL_H Oh ( r , x, y) := lin GEN_SOL_SEP f— , 

y-^y/x l x 


1 - k 

(1 - k) q, x, y ) 


( lin r) - y, x, y 
y-*x-y 


GEN_SOL_EX(p, q, x, y) : = J p dx + J |q - — J p dxj dy = c 
GEN_INT_FCT(p, p, q, x, y) := GEN_SOL_EX(p p, p q, x, y) 


J t du 

INT_FCT(t, u) := e 

GEN_FRAC_AUX (r , p, q, x, y) := lin GEN_SOL_HOM( lin 

[y^ x]-> [y - q, x - p] [y, x]-»[y + q, x + 

r, x, y) 
p] 
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#13: GEN_LI N_FRAC(r „ a, b„ c, p„ q„ k, x„ y) := GENJRAC_AUX r. 


c -p - a -k 


b k - c -q 
a -q - bp 


a -q - b p 


#14: AUX_LAGR(p, q, x, y) := 


I ./ 

J q e 


J p dx 

q - e dx + c 


#15: du 

LAGR(p, q, x, y , u, cl) := y = lin (xp + q), lim p = cl, x = AUX_LAGR . 

u->cl v-*cl p - y 

d d d ] 

— q — p — q 

du du du 

- , u, x , y = AUX_LAGR , - , u, x p + q 

p - u [ p - u p - u 

Example 7: Find the general solutions of the given differential equations 
and verify them. Find also the special solutions containing the 
given points . 


a) x-y' - 4y + 2x +4 = 0; 


P (1 I -3) 


b) y'-cos x + y-sin x = tan x; 


P (71/6 I 1) 


c) ( 1 — x ) y' =1 + x-y; 


P (2 | 2) 


These linear monic differential equations must be brought into the form/+ p(x).y = q(x). Then use 
LINEARl from ODE1 .MTH or GEN_LIN1 (p,q,x,y) from ODE 1_EXT.MTH to find the general 
solution written with a single symbolic constant c. 


Example 7a 


f a _ 2 a r\ t ^ 2x + 4 

x y ~4y + 2x +4 = 0 y'--y = 

x x 


2 x + 4 


#16: GEN_LIN1 


4 2 

#17: y = cx +x +1 

Check for correctness of the solution 

d 4 2 4 2 2 

#18: x (e x + x + 1) - 4 (c x + x + 1) + 2-x +4 = 0 

dx 

4 2 

#19: -3 = c 1 +1 +1 

4 2 

#20: S0LUE(-3 = c 1 +1 + 1, c. Real) 

#21: c = -5 


Special solution: y(x) = -5x 4 + x 2 + 1 


#22: LINEARl 


2 x + 4 4 2 

, x, y, 1, -3 = (y = - 5-x + x + 1) 


using LINEARl gives the special solution promptly. 
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„ - , . , sinx tanx , tanx 

Example 7b y cosx + jsmx = tanx <^> y+ y = <^> y+tmx-y = 

cosx cosx cosx 


#23 : 


#24: 


#25: 


LINEARl 


TAN(x), 


SIN(x) 


y = COS(x) 


COS(x) 

y0 


COS(x0) 2 

2 COS(x0) 


2 COS(x) 


GEN_LIN1 


TAN(x), 


SIN(x) 


#26: y = c-COS(x) + 


COS(x) 

1 

2 COS(x) 


#27: 1 = 


sT3 c sT3 


Which solution is more comfortable? #24 or #26? 


vT3 c ^ 

#28: S0LUE|1 = + , c. Real 


#29: c = 


2 W3 


2 _ ^ 2 

#30: y = I 1 COS(x) + 

3 3 J 2 COS(x) 


□r directly with LINEARl (...) 

LINEARl 


#31: 


SIN(x) n 

TAN(x), , x, y, , 1 

2 6 


Example 7c 


COS(x) 

(i-x 2 )/ = i+xj <^> /+ 


2 sf3 2 


-COS(x) + 

3 3 J 2 COS(x) 


' X ^ 


1-x 2 


y= r? 


#32: 


LINEARl 


#33: y = 


2 2 
x - 1 1 - x 

2 2 2 
LN(^(x - 1} + x> - LN(^(x0 - 1} + x0> - y0 ^(x0 - 1) 


sT(x - 1) 


#34: 


GEN_LIN1 


2 2 
x - 1 1 - x 


#35: y 


c - LN(^(x - 1} + x> 


^(x - 1} 


#36: y{x) := 


c - LN(\T(x - 1> + x) 


vffx - 1) 


For verfication of the result we define j;(x). 
For proceeding in this session we ’’redo” this 
definition and go on. 


#37: 

a - 

2 

x > y J (x> 

- a 

+ x y(x>) 

#38: 

0 




#39: 

(1 - 

2 

x > y J (x> 

= i 

+ x-y(x) 

#40: 

SF<X> 

: = 



#41: 

i, : = 
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In later Derive versions LINEAR1_GEN() was included in 0DE1.MTH together with other 
_GEN()-functions in order to create general solutions using a symbolic constant (c by default). 


#42: LINEAR1_GEN 


X 1 


2 

c - LN(\T{x - 1) + x) 

2 

2 


9 - 

2 

L x - 1 

1 - x , 


^(x - 1) 


and finally the requested special solution 

X 1 


#43: 


LINEAR1 


2 

x - 1 


x, y, 2, 2 


1 - x 


#44: y = 


2 - sT3 - LN{{2 - 13) <I(x - i) + x)) 


sT(x - 1} 


Both of the following equations are of the Bernoulli-type: 

y' + p(x) • y = q(x) • y k (k = const.) 


In ODEl . MTH we found BERNOULLI (p,q,k,x,y,xO,yO). GEN_BERN (p, q, k) gives the 
general solution with constant c. In DERIVE5' s ODEl . MTH from today we can find 
BERNOULLI ODE () and BERNOULLI ODE GEN ( ) . 


Example 8 : 


x 2 

a) y'-y = e - y ; 

P (0 | 2) 

3 2 2 3 

b) y'-x + y-x - (x +1) -y ; 

verify the solution 


BERNOULLI() now gives the 
BERNOULLI number! 


#45: BERNOULLI (1, e , -1} = 


In 1991 Bernoulli ( ) delivered the solution of the differential 
equation: 


1: BERNOULLI (1, e , -1) 


2 2 e - 2 x 

2: y = + e 

3 


2 2 xO 

yO e 


3 xO n 


2 e 


Example 8 a 

This is DERIVE5 and ODEl EXT. MTH 


COMMAND: 

Enter op- 
User 


rjFtajTSn Build Calculus Declare Expand Fact 
Options Plot Quit Remove Simplify Transfe 
i on 

Free : 10 


#46: BERNOULLI _ODE( 1 , e , -1) 

x 


2 2 e - 2 -x 

#47: y = + e 


2 a 2 x0 
y0 e 


2 e 


3 x0 


#48: GEN_BERN(1, e , -1> 

x 

2 2 e - 2 x 

#49: y = + ce 
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#50: BERNOULLI _0DE_GEN<1, e , -1) = 


2 2 e - 2x 

y = + c e 


#51: BERNOULLI _0DE<1, e , -1, x, y, 0, 2) 

0 

2 2 e - 2 0 

#52: 2 = + c e 


x - 2 -x 

2 2 e 10 e 


#53: SOLUE 


2 2 e a - 2 0 

2 = + c e , c. Real 


#54: c = 


10 


Compare with the Derive solution from above! 
Example 8b 


#55: 


#56: 


GEN_BERN 



2 

x + 1 


3 


x 


4 2 

1 2cx +2x + 


2 2 
y 2 -x 


3 


1 


#57: 


#58: 


1 


2 

y 


y(x) := 


1 

4 2 

2-c-x + 2 x +1 

2 

2 x 


>T 


2 

2 x 


4 2 

2-c-x + 2 -x 


+ 1 


3 2 2 3 


#59: 


(x)-x + y(x)-x - (x + l)-y(x) = 0 

#60: 

0 

= 0 

#61: 

y 

: = 

#62: 

y 

= y 


Example 9: Test the following differential equations on homogeneity. 

Find two ways for obtaining the general solutions . 

. y + x 

a) y =- ; P (1 1 2) 

y-x 

v e y/x — x 

b) y T = — — ; P (1/2 | — 3/2) Verify the solution! 


Derive 2.05 had a function homogeneous_test ( ) , which disappeared in later versions. 

The equation is homogeneous /7homogeneous_test ( ) gives back an expression which is free of 
the independent variable. We can now reanimate this function or make a new one, which uses the 
definition of homogeneous expressions: substituting k x and ky for x and y should result in an equiva- 
lent expression. 
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H0M0GENE0US_TEST (r, x, y) := lim r 

y-»x v 


H 0M0GENE0U S _T ES T 


H 0M0GENE0U S _T ES T 


y + x 


V - x 


y + 1 

y - i 


^y/x 

ye - x 


y/x 


y - e 


H 0M0GENE0U S _T ES T 


y/x 

ye - 1 


y/x 


= y 


my_hon_test (r, x, y) := 
If r - SUBST (r, [x, y], 

''homogeneous'' 

''non homogeneous'' 
''non homogeneous'' 

. [k x, k y]) = 0 

my_hom_test | 

[ y + x 1 


i — Jiumuyciicuus 

l y - x J 

my_hom_test 

r y/x 

ye - x 

= homogeneous 

y/x 

L x e 

my_hom_test 

r y/x 

ye - 1 

= non homogeneous 

y/x 

L x e 


Proceed with solving the equation: 

y + x 


#63: HOMOGENEOUS 


y - x 


#64: 2 

LN{- x0 - y0 (2 x0 - y0)> 


LN 


2 2 2 1 

x0-{x +2x-y-y) 


LN 


X0 


This was the 1991 -recipe for manipulating this expression in order to find the solution: 

Build #64, ENTER, *, 2, ENTER, E (xp) , R(ecip), A ENTER 

Manage Substitute xO =1, yO = 2, A ENTER 

2 2 

x +2 x y - y =1 

Let's do the same without the "Building procedure" from Derive-for-DOS-times: 

#64 * 2 


#65: 


LN 


LN(- x0 - y0 (2x0 - y0)> 


2 2 2 
x0 (x + 2- x- y- y) 


#66: LN{- x0 - y0 (2x0 - y0>) - LN 


2 2 2 
x0 (x + 2 - x - y — y } 


LN 


X0 


2 LN 


x0 


2 
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e n (#66) 


2 2 2 2 2 
#67: ^LN(- x0 - y0 (2 x0 - y0)> - LN(- x0 (x + 2 x y - y )/x ) = - 2LN(x0/x) 

e 


# 68 : 


2 2 2 
x (x0 + 2-x0-y0 - y0 ) 

2 2 2 
x0 (x + 2 - x - y - y ) 


X0 


1/#68 


2 2 2 
#69: x (x0 + 2x0y0 - y0 ) 


2 2 2 
x0 -(x + 2 - x - y - y } 

2 2 2 
x0 (x + 2 x y - y ) 

2 2 2 
x <x0 + 2x0-y0 - y0 ) 


x0 

2 


#70: 


x0 


#70/(x0 A 2/x ft 2) 


#71: 


2 2 2 2 
x0 (x + 2-x-y - y ) x0 

2 2 2 2 
x (x0 + 2x0-y0 - y0 ) x 


x0 


#72: 


2 2 
x + 2 -x -y - y 

2 2 
x0 + 2 - x0 ■ y0 - y0 


Substitute forxO = 1 , yO = 2, leading to 

2 2 
#73 : x + 2-x-y - y =1 


#74: GEN_SOL_HOM | — — — , x, y 

y - x 


#75: 


#76: 


LN{- x - y (2 x - y)) 


#77: - 


LN{- x - y (2 x - y» 


2 


#78: {- LN{- x - y (2 x - y))/2 = - n i/2) (-2) 


#79: - x - y (2 x - y) = -1 


If you exactly know what to do, then it is easy in Derive5 , too: 

2 2 

SUBST ( l/e A (2 *#64) / (xO A 2/x A 2) , [xO, yO] , [1,2]) which results in x +2-x-y-y =1 


#80 : H OMOGENEOU S _GEN 


(t^ 


#81 : - LN 


sT(- x - y (2 x - y>) 


= LN(x) + c 


#82: HOMOGENEOUS |— — „ x, y, 1, 2 

y “ x 


#83: 


LN 


2 2 
x + 2-x-y - y 


= LN(x) 


Transform to find a nice reprecentation of the solution! 
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#84: HOMOGENEOUS 


V/x 

ye - x 1 3 

J- X, y j- r 

y/x 2 2 


-3 y/x 

#85: e - e = LN(2 x) 

-3 y/x 

#86: SOLUE(e - e = LN(2 x), y, Real) 


#87: y = x LN(1 - e LN(2 x>) - 3 x 

Or do it stepwise: 


#88 : H OMOGENEOU S _GEN 


y/x 


^y/x 

ye - x 


y/x 


#89: 

- e = LN(x) + c 


#90: 

f f 

SOLUE SUBST - e = LN(x) + c. 

[x, y]. 

#91: 

-3 

c = LN{2) - e 


#92: 

y/x 

S0LUE{- e = LN{x) + LN(2) - e 

-3 

, 9> 

#93: 

3 

y = x LN{1 - e LN{2 x>) - 3 x - 

2 nix 


v- - -f ])• •] 


V y = X LN(1 - e -LN{2 x)) - 3 x 

As we didn't ask for the real solutions only, we have to select the real one: 

3 

#94: y = x LN{1 - e LN{2 x>) - 3 x 


Verification of the solution: 


#95: y := x LN(1 - e LN(2 x)) - 3 x 

^y/x 

d y -e - x 

#96: — y = 0 

dx y/x 

x e 


The counter example: 


H OMOGENEOU S_GEN 


y/x 

ye - 1 


y/x 


inapplicable 


Example 10: Compare the integral curves of 

2 2 0 

y -x 2xy 

y f = and y , = — ? y intersecting in (2 | — 1 ) 

2xy y -x 

Plot both curves together with their tangents in P 
on the same axes. 
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#99 : GEN_SOL_HOM 


2 2 1 
y - x 


2 x y 


2 2 
x + y 


#100: - LN 


2 2 

#101: (- LN((x + y )/x) = c) 

e 

2 2 

x + y — c 

#102: = e 


#103: HOMOGENEOUS 


#104: - LN 


2 2 
V - x 


2 x y 


x, y, 2, -1 


2 2 
4 <x + y ) 

- LN f— — — 1 

2 

- LN l 2 ] 

5 -x 



#105: SOLUE 


# 106 : y 


- LN 


2 2 
4 (x + y > 


5 -x 

vT2 vTx vT(5 - 2 x) 


LN| — |, y 


vI2 vTx vT{5 - 2 x) 


This turns out to be a family of circles with centers on the x-axis. The solution containing point P is 
one of the semicircles from #106. The direction field for only one single point gives the tangent of the 
integral in this point. (Plot #106 - #109 and see the circle, the point and the tangent (-2 <t< 2). 


#107: FELD 


2 2 
y - x 


2 xy 


x, y, 2, 2, 1, -1, -1, 1 



4 t 3 t 

+ 2 , 1 


#111: 


, x, y, 2, -1 


2 2 

LN(y (y - 3 x )) LN(ll) 

#112: + LN(x) + LN{2) = LN(x> - LN(2) 

3 3 

2 2 

#113: {- LN(y (y - 3 x ))/3 + LN(x) + LN(ll)/3 - LN{2) = LN(x) - LN(2>) 3 

e 


#114: 


#115: 


#116: 


11 x 


2 2 
8 y (y - 3 x ) 


11 x 


2 2 
8 y (y - 3 x ) 


2 2 
8 y (y - 3 x ) 


2 2 

ll = y (y - 3 x ) 


In 1991 we couldn't do implicit plots. Now we can. In DNL#3 we solved equation #116 for y and re- 
ceived three solutions (= three explicit branches of the integral cuve). This is the plot from DNL#1. 
There is a third part which is outside of this plot window. 
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Without any doubt I believed 
Derive, but today I know a bit 
more: 

Check if both points (2\-1) and 
(-2\1) really fulfill equation #1 16? 


The final plot shows both curves together with the tan- 
gents. And we see that Derive 5 delivers another plot. 
Solving #116 now for y and plotting the three branches 
gives exact the same figure - and checking the points 
works, too. 



Example 11: Find the curves which intersect the circles x + y = r at 
angles of tu/4 . 

Plot some circles and a family of curves representing the 
general solution of the underlying differential equation. 


This problem leads to the homogeneous DE y r = 


y-x 

y + x 


We start working with rectangular coordinates and switch later to polar coordinates. 
Compare my "do it yourself’ GEN_SOL_HOM with professional HOMEGENEOUS_GEN: 

#122: GEN_SOL_HOM ' ^ 


y + x 


#123: 


ATAN 


LN(sT(x 


2 

y )) 


#124: H OMGGENEOU S _GEN 


y - x 


y + x 


LN 


^(x 


2 

y ) 


#125: 


- ATAN 


= LN(x> + c 



These partially plotted spirals are not really satisfying me. We switch to polar coordinates substituting 
for x = r cos(cp) and for y = r sin(cp). 

#127: r :e Real (0, ®) 

2 2 

[ r SIN(ip) ) LN((r COS(ip)} + (r SIN(ip)} > 

= c 

r COS(ip) J 2 


#129: 

#130: 

#131: 

#132: 


f -P 1 1 

n FLOOR + - LN(r) - ip = c 

In 2 J 

f 11 

n FLOOR + - LN(r) - ip = - LN(c> 

l n 2 J 


SOLUe|h FLOOR + -^-J 


- LN(r) - ip = - LN(c), r 


-FLOOR(ip/n + 1/2) - ip 

v = c -e 


#134: - LN(r) - ip = c 
#135: - LN(r> - ip = - LN(c) 

#136: SOLUE(- LN(r) - ip = - LN(c), r) 

“0 

#137: r = c e 


#138: UECTOR(c e c, 0-5, 5, 0-5) 
#139: UECTOR(c, c, 0.5, 5, 0.5) 
#140: 2 


-o 

#141: 2 e 
#142: [2, 0] 


I y - x 

, x, y, 2, 2, 1, 0, 0, 1 

y + x 


#144: 


^2 t 

+ 2 , 

2 


sf2 t 
2 


#145: [2, t] 


p 26 


Josef Bohm: Solving ODEs using DERIVE 


D-N-L # 3 



A family of circles together with a family of integral curves. The tangents in 
point (2|0) are plotted 


Example 12 : A point moves on a curve in x-y-plane in such a way that the 
angle formed by the tangent of the curve and the x-axis is 
three times the angle between the radius vector and the x-axis. 
Find the Cartesian equation of the family of curves satisfying 
this condition. 

Plot a family of curves, and plot the special curve containing 
point P (-3 | -2 ) . 


Using / = tan(a) = tan(3B) and tan(G) = y/x you find after some calculations - let Derive do this job - 


applying trigonometric identities the homogeneous differential equation y' = ~ . 

x -3xy 


#185: TAN<3 p) 

#186: Trigonometry := Expand 


#187: 


#188: 


4 SIN ATAN 


4 SIN(p) COS(p) 
2 

1 - 4 SIN(p) 
COS I ATAN 


SIN(p) 


i)F ■“(•“[t-]) 


COS(p) (4 SIN(p) 
SIN ATAN 


1 > 


d) 


4 SIN ATAN 


-IF 


COS ATAN 


2 2 

y (3 x - y ) 


)] 


4 SIN ATAN 


IF 


# 189 : 


2 2 
x (x - 3 y > 
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This is the original file from 1 991 : 


#16: GEN_SOL_HOM 


2 3 

3 x y - y 


3 2 

x - 3 x -y 


#17: - LN 


2 2 
x + y 


#18: - LN 


\Tx ■ \Ty 

2 2 
x + y 


\Tx ■ \Ty 


- LN(c) 


2 2 
x + y 

#19: = c 

\Tx ■ \Ty 


# 20 : 

#21: 

# 22 : 

#23: 

#24: 

#25: 


lots? 


o we use Polar coordinates! 


2 2 
(r COS(ip)) + (r SIN(ip)) 

sT(r COS(ip) W(r SIN(ip)) 

r 


vT{COS(ip)W(SIN(ip)) 

2 


= c 


COS(tp) -SIN(tp) 
2 


2 r 


SIN(2 ip) 


= c 


This is Bernoulli's Leminscate! 


#26: r 


vf2 c vT(SIN(2 ip)) 
2 


#27: 


#28: 


#29: 


#30: 


y[2 c vT(SIN(2 ip)) 

r = 

2 

2 2 
(-3) + (-2) = 

^(-3) ^(-2) 

13 - \T6 

- = c 

6 

13 - \T3 - \T(SIN{2 ip)) 

r = 

6 


The following commands are the base for the plot. 

It is important to switch between rectangular (#31, #32) 
and polar coordinates (#33). 

It is one of Derive 's unique features that one is able to 
have both coordinates systems available in the same plot. 


For plotting the family of lemniscates set -n/2 < cp < nil. 


#31: FELD 


2 3 

3 x y - y 


3 2 

x - 3 -x y 


x, y, -3, -3, 1, -2, -2, 1 


#32: 


9 - \T13 - 1 46 - \T13 - 1 

3, 2 


169 


169 


^ c ^{SIN{2 ip)) 

#33: U ECTOR | , c, 0.5, 5, 0.5 


# 1 : 


#2: 


HOMOGENEOUS 


2 

3 -x y 


3 


3 

y 

2 


LN 


x - 3 x y 

2 2 2 1 
54 (x + y ) 

3 

169 x y 


2 


x, y, -3 



-2 


n - i 


e*(#2): 

13 sT6 


2 2 1 

#3: 18 (x + y W 

3 

l * -y . 


X 

3 


This cannot be plotted in Derive 5, but interes- 
tingly enough it can be plotted in Derive 6. 

And here we find one part of the lemniscate. 



#4: 


[-3. -2] 
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Finally we plot a family of the complete lemniscates together with the tangent in P(-3|-2) and the ray 
from 0 to P. We will show that the relation between the angles given in the task is true. 



Example 13: Show that the given differential equations are exact, give 

their general solutions and if there is a point given find also 
the respective special solution. 

a) 2x 2 - y 2 + y - y' (2xy - x + 4y) = 0; P(l/-1.5) 

b) 2x + e X in y + (e X y')/y = 0 

c) 2xy-y'(2x 2 +y)=0; Q ( — 2 | 3 ) 


DEs of form p(x,y) + q(x,y) / = 0 may be exact. Exactness seldom is obvious, so 

EXACT_TEST (p, q, x, y ) or EXACT_IF_0 (p, q, x, y) in earlier versions will test this attribute. 


EXACT_TEST is not included in recent ODE1.MTH. The test has been included into the EXACT()- 
f unctions for solving the DEs. For didactical reasons it might be useful to ask students to produce a 
selfmade EXACT_TEST-function. In the following I give the original EXACT_TEST from 1991 and an 
updated form for DERIVE 5 (and 6, of course). 

d d 

#246: EXACT _TEST (p, q, x, y) := — p q 

dy dx 

EXACT _TEST_N(p„ q, x, y) := 

If DIF(p, y) - DIF(q, x) = 0 
#247: "exact" 

"not exact" 

"not exact" 


2 2 

EXACT _TEST(2 x - y + y, - (2 x y - x + 4 y)> = 0 
2 2 

EXACT _TEST_N(2 x - y + y, - (2 x y - x + 4 y)) = exact 


EXACT _TEST 


2 x 


+ e X LN{y), 



= 0 


EXACT _TEST_N 


2 x 


+ e X LN{y), 



exact 


2 

EXACT_TEST(2 x y, - 2 x - y) = 0 

2 

EXACT_TEST_N{2 -x -y, - 2 x - y) = not exact 
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2 2 

#254: EXACT (2 x - y + y, - (2 x y - x + 4 y)) 

3 2 3 2 

2x 6-y + 2 x0 + 3-x0-y0-(l - y0> - 6y0 

#255: + x y (1 - y> 


= 0 


2 2 3 

#256: EXACT |2x - y + y, - (2 x y - x + 4 y), x, y, 1, — 

3 2 

2 x 24 y - 91 

#257: + x y (1 - y) = 0 

3 12 

3 2 

#258: 8 x + 12 x y (1 - y> - 24 y + 91 = 0 

2 2 

#259: EXACT_GEN(2 x - y + y, - (2 x y - x + 4 y)) 


2-x 2 

#260: + x y (1 - y) - 2 y = c 


#261: 


91 

12 


leads to the same solution (#258) 
Example 13 b 


#262: EXACT _GEN 


2 x + 0 X -LN(y), — 


x 2 

#263: £ LN(y) + x = c 

x 2 

#264: SOLUE(e LN(y) + x = c, y) 


#265: 


y = e 


#266: SOLUE 


-x 

2 

- (c 

- x ) 


r X 1 



x e 


EXACT 

2 x + e LN(y>, 

i- 9 


L 9 , 



— x x0 2 2 

#267: I x0 - x e (e LN(y0> - x + x0 > 

y = I F [e 1 1 v y0 > 0, e 


Example 13 c 

This equation is obviously not exact (see above). We assume, that we know an integrating factor 
p = 1/y 3 . Now we are able to use USE_INTEG_FCTR (p, p, q, x, y, xO , yO ) . ( This function is not 
available in later Derive versions.) 
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#268: USE_INTEG_FCTR(p, P, q, x, y, x0, y0> := EXACT(p p, p q, x, y, x0, y0) 
USE_I NTEG_FCTR 


#269: 


1 2 
, 2x-y„ -2-x - y 


2 2 2 2 
y0 -x - y{y{x0 + y0> - y0 ) 

#270: = 0 

2 2 
y0 y 


#271: 


USE_INTEG_FCTR 


1 2 

, 2 xy „ - 2 x - y, x, y, -2, 3 


9 x - y (7 y - 9) 

#272: = 0 


#273: 


9 y 
GEN_INT_FCT 


1 2 
, 2-x-y, -2 x -y 


x 1 

#274: + = c 


2 y 


#275: = c 

9 


o , • X 2 1 7 

Solution: —r H — = — 

/ 7 9 


In times of Derive 5 and higher all the auxiliary functions (see later in example 14) are implemented in 
on function to solve immediately differential equations, which are supposed to be soluble by applying 
an integrating factor: I recommend to compare the recent ODE1.MTH file with the auxiliary procedures 
which where necessary in 1991. 1 can imagine that it might be useful for students studying differential 
equations to follow the old "recipes" (which can support understanding the solving process). 


#277: INTEGRATING^ ACT OR (2 x y. 


9 -x - y (7-y - 9) 

#278: = 0 


9 y 


2 

2-x - y, x, y, -2, 3) 

Compare with #275 


File ODE1_1991.dfw is the original ODE1.MTH from Derive 2.05 together with 
ODE1 EXT.MTH from DNL#3. 


Example 14: Try to find a solution of the DEs given below using the auxil 
iary functions of ODE1.MTH and ODEl_EXT.MTH concerning the 
application of integrating factors. Give the integrating fac 
tors . 

a) 2xy - y' (2x 2 + y) = 0 

2 2 4 

b) (x + y)(xdy-ydx) = (a + x) x dx 

c) (x cos y - y sin y) dy + (x sin y + y cos y) dx = 0; P ( 1 | tt/2 ) 

d) (2x 3 y 2 - y)dx = (x - 2x 2 y 3 )dy 
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The DEs must be brought into the form p(x,y ) + q(x,y) / = 0. 

DERIVE can be helpful to find an integrating factor p. 

Preloaded is File ODE1 .MTH, copyright (c) 1 990 by Soft Warehouse, Inc. 

& ODE1_EXT from Derive Newsletter #3 

2 3 

#41: FREE_OF_I NDEPENDENT_TEST ( 2 - x - y „ - 2 x - y) = 

9 

simplifies to an expression free of x. Hence we get the solution by means of the next function: 


#42: FREE_OF_I NDEPENDENT (2 - x - y , - 2 x - y) 

2 2 2 2 
y0 x - y(y(x0 + y0) - y0 ) 

#43: = 0 

2 2 
y0 y 

INT_FCT f — „ y| = — 

#44: ( y J 3 


Follow example 1 3c from above 


Example 14 b 

The given equation must be rewritten as 


{x +y ){xy' -y) = (a + x)x 
-x 2 y-y 3 - ax 4 -x 5 +(x 3 +xj^ 2 )^' = 0 
x 2 y + y 3 + ax 4 +x 5 +(-x 3 ~xy 2 )y r = 0 


5 4 2 3 3 2 

#45: FREE_OF_I NDEPENDENT _TEST (x +ax + x -y + y , - x - x-y } 


#46: - 


2 2 
4 (x + y ) 


5 4 2 3 

x + ax + x-y+y 


which is not free of independent x, but 

5 4 2 3 3 2 

#47: FREE_OF_DEPENDENT_TEST(x +ax + x -y + y , - x - x-y ) 


#48: - 


which is free of dependent y. This leads to 

5 4 2 3 3 2 

#49: FREE_OF_DEPENDENT (x +ax +xy+y,-x -xy) 


2 3 

x 9 9 

#50: + a x - - 

2 x 3 

3 -x 


5 2 3 

3x0 - 6 -x0 -y0 - 2-y0 

a x0 - = 0 

3 

6 x0 


INT_FCT I — , xl = — — 

#51: l x 1 4 


intesratins factor!! 


#52: 


GEN_I NT_FCT 


1 5 4 2 3 3 2 

, x +a-x +x-y+y,-x -x-y 


2 3 

x y y 

#53: + a x - - = c 

2 x 3 

3 -x 
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Now in DERIVE5: 


5 4 2 3 3 2 

#279: I NTEGRAT I NG_FACTOR{x + a x + x y + y , - x - x y ) 


2 3 

x 9 9 

#280: + a x 

2 x 3 

3 -x 


5 2 3 

3-x0 - 6 x0 -y0 - 2-y0 

a - x0 - = 0 

3 

6 x0 


Example 14 c 

After rewriting we proceed: 


#54: FREE_OF_DEPENDENT_TEST(x SIN(y) + 9 COS(y), xCOS(y) - 9 SIN(y)) 

#55: 1 

#56: FREE_OF_DEPENDENT |x SIN(y) + 9 COS(y), x COS(y) - y SIN(y), x, y, 1, 

x 

#57: e (y COS(y) + (x - 1) SIN(y)} = 0 

#58: INT_FCT (1, x) = e* 

#59: GEN_INT_FCT(e X , x SIN(y) + y COS(y), x COS(y) - ySIN(y)) 

#60: e (y -COS(y) + (x - l) -SIN(y)> = c substitute for X and V 


Please check the result with the respective Derive5 function! 


Example 14 d 

2 2 

32 23 4 x (y - x ) 

#62: FREE_OF_INDEPENDENT_TEST<2 x y - y, 2 x y - x> = 

3 

2 x y - 1 
2 2 

32 23 4 y (x - y ) 

#63: FREE_OF_DEPENDENT_TEST(2 x y - y, 2-x -y - x) = 

3 

2 x y - 1 


neither free of x nor of y ; we test other possibilities for the integrating factor: 


We proceed in DNL#3 (1991) - style: 

32 23 2222 

#64: EXACT_TEST(x y (2 x y - y), x y (2 x y - x)) = 6 x y (x - y > 

xy is no integrating factor 
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3 2 2 3 4 3 2 

#65: EXACT _TEST{{x + y) (2 x y - y), (x + y) (2 x y - x)) = 4 x y + 6 x y - 

2 3 4 

6 x -y + x (1 - 4y ) - y 

no luck with (x + y) 

22 32 22 23 2222 

#66: EXACT_TEST(x y (2 x y - y), x y (2 x y - x)> = 6 x y (x - y ) 


no success with x*2 y*2 

l 


#67: 


EXACT _TEST 


2 2 
x y 


3 2 1 2 3 

(2 x y - y), (2 x y - x) 


2 2 
x y 


= 0 


This was the right guess! Let's finish the solution 

GEN_INT_FCT 


# 68 : 


1 3 2 2 3 

, 2x -y - y, 2-x -y - x 


2 2 
x y 


3 3 

x -y + x y +1 

#69: = c 

x y 


with MONOMIAL_TEST (p, q, x, y) DERIVE can find an integrating factor of the form x m y n if there 
exists one. 

This is the way we do it now: 

3 2 2 3 

#281: INTEGRATING_FACT0R_GEN(2 x y - y, 2 -x -y - x) 

#282: inapplicable 


3 2 2 3 

MONOMI AL_TEST(2 x y - y, 2 x y - x) 

#283: 


#284: 


I NTEGRAT I NG_FA CT OR_GEN 


1 3 2 

(2 x y 

2 2 
l k y 


2 2 
x y 

1 2 3 

y), (2-x y - x) 

2 2 
x y 


3 3 

x -y + x -y + 1 

#285: 

x y 


Example 15: (x-2y+5)dx+ (2x-y+ 4)dy = 0 

a) Find the general solution using both ( x 0 | y 0 ) and constant c. 

b) Find the special solution with x = 1 , y = 1 . 

c) Sketch the integral curve and the direction field. 


LIN_FRAC ( r , a, b, d, p, q, k, x, y r xO, yO) simplifies to an implicit solution of a 
linear fractional equation y' = r-((a-x+b-^+d)/(p-x+q- 7 +k)). However if q-a - p-b = 0, instead use 
FUN_L I N_C C F described above. (Online - Help). 

We will tackle this problem in the Derive 5/6 mode only and add the 1991 plot for c) 


f -x + 2-y - 5 

#286: LIN_FRAC , -1, 2, -5, 2, -1, 4 

( 2 x - y + 4 


LN 

x + 1 


x - y + 3 

- LN 

sT 

r x0 - y0 + 3 V 




#287: 

X0 + 1 


3 



3 




(x + y - 1) J 



L (x0 + y0 - 1) 


= LN(x + 1) - LN(x0 


+ 1 ) 
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#288 : LIN_FRAC_GEN 


-x + 2 - y - 5 
2 x - y + 4 


-1, 2, -5, 2, -1, 4 


# 289 : 


LN 


(x - 1) (x - y + 3) 


<x + y - 1) 


= LN(x + 1) + c 


#290: 


I x + 1 1 - %r 


x - y + 3 


#291: LIN_FRAC 


<x + y - 1) 
-x + 2 - y - 5 


= e (x + 1) 


2 -x - y +4 


-1, 2, -5, 2, -1, 4, x, y, 1, 1 


#292: 


LN 


(x - 1) (x - y + 3) 


12 (x + y - 1) 


= LN 


x + 1 


sT3 - |x + 1 1 - %T 


#293: 


x - y + 3 


(x + y - 1) 


SIGN((x + y - 1) (x - y + 3)) 


x + 1 
2 



This is the graph of #293. Building e A (2*#292) 
leads to expression #294 which can be plotted 
together with the direction field. As you can 
see a discontinuity seems to appear at (-1\2). 

(Let students deduce this point from the task!) 

In 1991 we received the same solution - for- 
tunately - but we couldn 't do the implicit plot. 
So I solved #294 for y and plotted the soluti- 
ons together with the results of my FELD- 
function. 


x - y + 3 

= 1 

#294: 3 

3 (x + y - 1) 


f -x + 2-y - 5 

#295: DIRECTIONS I ELD , x, -3, 3, 24, y, -2, 3, 20 

( 2 x - y + 4 


The right graph is "Derive for DOS made". 






We bring the equation into the form p(x v-y) = q(y) and then use CLAIRAUT (p, q, x, y , v, cl ) 
taking v as a variable instead of/. 


#296: 

#297: 

#298: 

#299: 

#300: 

#301: 


#302: 

#303: 
#304: 
#30G : 
#306: 
#307: 

#308: 

#309: 

#310: 


2 

y = xu + u - u 
2 

x - u — y = u - u 

2 

CLAIRAUT (x -u - y„ u - u) 

[ 2 1 

Lcx-y-c +c=0, x-2u+l=0J 

f x + 1 

SOLUE(x - 2-u + 1 = 0, u> = u = 

f 2 x + 1 

SUBST y = xu + u - u , u, 

l 2 

2 

(x + 1) 

9 = 

4 

2 

SOLUE(cx -y-c + c = 0, y) 
y = c(x - c +1) 

UECTOR(y = c (x - c + 1>, c, -3, 3, 0.5} 
[2, 1] 

SOLUTIONS (GUEST (y = c (x - c + 1>, [x, y], 

^ 3 3 sf5 

2 2 ' 2 2 


In order to find the singular solution 
we solve the 2 nd component of #299 
for v and then substitute for v in the 
given differential equation (#296). 


This parabola is the singular solu- 
tion. 

The 1 st component of #299 gives 
the general solution. 


[2, 1]>, c) 


U ECTOR y = c (x - c + 1}„ c. 



sfS 

2 


]) 


9 


(^ + 3) (2 x - >15 - 1) 

4 


9 


(3 - ^IS) <2 x + >15 - 1) 
4 



The special solution(s) are the two 
tangents from point P to the singu- 
lar solution. 
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Example 17: Find a curve with its tangent's intercept between the axes 
having constant length a = 2 . 

Give both general and singular solutions. 


The problem leads to the differential equation y = xy'± . y , which is of Clairaut-form. 


V( 1 + t 2 ) 


#311: 


cl : = x - u — y = 


2 u 


#312: 


#313: 


CLAIRAUT 


sl<l + u ) 
2 u 


x u - y. 


sT{l + u ) 


2 2 

c-x-\T(c + 1) - y-\T(c + 1) - 2 c 


#314: 


SOLUTIONS 


2 

s[(c +1) 

2 

2 3/2 

{v + 1) 


= 0, x 


2 3/2 

+ i> 


= 0 


= 0, U 


#315: 



2/3 2/3 


2/3 2/3 



2 - X 

, - ^ 

2 - x 


2/3 

2/3 





L X 


X 



#316: U ECTOR 




2/3 2/3 


2/3 2/3 



cl, u. 


2 - x 

, - %T 

2 - x 



2/3 

2/3 








x 


X 




#317: 



2/3 2/3 


2/3 2/3 


2/3 2/3 

x - \T 

2 - x 

2/3 1/3 

- y = 2 x - 

2 - x 

"7 

X 

1 

■■ 

2 - x 

2/3 

2/3 

2/3 





x 


X 


L X 


2/3 1/3 

2 x ->T 


2/3 2/3 

2 - x 


2/3 


2 2/3 2/3 3 

#318: y = (2 - x ) 


Derive 2 was unable to solve #5 for v (see below), so I tried to find another way: 


2 3/2 

x (u +1) - 2 

#5: = 0 

2 3/2 

+ 1) 

Substitute for v in the second component and in the equation: v = tan(cp). 
Solve for x and y to obtain a parameter representation of the singular solution. 

2 

#6: x - 2 COS(ip) |COS((p)| = 0 

2 

#7: x = 2 -COS(ip) |COS((p}| 

2 2 TAN(m>} 

2 COS{(|>) |COS((p)| TAN((|») - y = 

#8: 2 

sT(l + TAN<(p) ) 

3 

#9: y = - 2 -SIN(ip) SIGN(COS{(p)) 

#10: [2 COS(ip) 2 |COS(ip)|„ - 2 SIN(*p) 3 SIGN(COS(m»>>] 


y = - 
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#318 ( DERIVE5 ) and #10 ( DERIVE 2) give an astroid and the gene- 
ral solution is a family of lines (tangents of the astroid). 


#319: SOLUE 


2 2 

e x -^(c + 1) - y-\T(c + 1) - 2 c 


= 0, y 


sT(c + 1) 


#320: y = 


c ■ (x ■ \T(c + 1) - 2) 


T(c +1) 


#321: U ECTOR 


c ■ (x ■ \T(c + 1) - 2) 

y = , c, -4, 4, 0.5 

2 

^(c + 1) 



Example 18: The following equations of form y = x p(y’) + q(y') are called 
LAGRANGE differential equations. 

Use L AGR (p,q,x,y,v,c) to solve the following equations: 

a) y = x (1 + y') + y' 2 

b) y' 3 - 3y' = y - x 

2 2 

c) (x y f + y) = y y' 

LAGR (p,q,x,y,v,c) with v instead of / simplifies to a vector of four components with the last two 
of them giving the parameter form of the general solution. If one is able to solve the 2 nd component for 
cl one will obtain the singular solution by substituting for cl in the 1 st component. LAGR ( ) is not part 
of ODE1.MTH but of ODE1EXT.MTH. Look at the examples: 

Example 18a 

r 2 -u -U 2 

#323: [y = x-(cl + 1) + cl , cl + 1 = cl, x = ce - 2 ■ (u - 1), y = ce -(u + 1) - u + 


No singular solution because of cl + l = cl. 

Integral curves + direction field: 

r 2 l 

#324: [ c -e - 2 (u - 1), c e (u + 1) - v + 2\ 


#325: UECTOR{ 


f 2 1 

fee -2-(u-l), c-e (u+l)-u +2j, c, -5, 5, 0. 


5) 


#326: SOLUE(y = x - (1 + u> + u , u> = 


\T(x - 4 x + 4-y) + x 


\T(x - 4 x + 4 y) - x 
u = v u = 


#327: DIRECTIONS I ELD 


#328: DIRECT IONS I ELD 


\T(x - 4 x + 4 y) - x 


x, -5, 5, 10, y, -4, 4, 8 


\T(x - 4 x + 4 y) + x 


-, x, -5, 5, 10, y, -4, 4, 8 
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Example 18b 


#329: LAGR(1, u - 3 u) 


#330: 


3 3 u (u +2) v (2 u + 3) 

y = x + cl - 3 -cl, 1 = cl, x = c + , y = c + 


#331: SUBST (y = x + cl - 3 cl, cl, 1) = (y = x - 2) 


singular solution! ! 


#332: U ECTOR 


3 u (u +2) u (2 u + 3) 

c + , c + 


, c , —3 , 3 , 0.5 


Example 18c 


2 2 f u x 

#333: SOLUE((x u + y) = y u, y) = y = v y 


sTu - 1 


sfy + 1 


Let's take the first solution: 


#334: LAGR 


#335: 


Jv - 1 
cl x 


y = 


#336: x 


^cl - 1 \Tcl - 1 

c 'V- 1 


cl c (sjy - 1) 

= cl, x = , y = c - \Tu 


#337: SOLUE 


(v-l I 




Example 18c 

After eliminating parameter v we find a family 
of hyperbolas as general solution. 

We have two singular solutions (lines). 

The result of the LAGRANGE DE derived 
from the second solution fory in #333 can be 
found on the next page. 


#338: y = 


#339: U ECTOR 


-, c , —3 , 3 , 0-5 


#340: SOLUTIONS 


c - x 
cl 

\Tcl - 1 
f cl -x 

#341: U ECTOR y = 

[ \Tcl - 1 

#342: [y = 0, y = 4 x] 


cl, cl = [0, 4] 


cl, [0, 4] 


#343: LAGR 


#344: 


sJu + 1 
cl -X 


0 


cl 


sTcl +1 sTcl + 1 


cl, X = 


- 2/sTu 

c-e -(sTu + 1) 2/sTu 

, y = c sTu e 
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